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ABSTRACT 

We present the results of an infrared (IR) study of high-redshift galaxy clusters with the MIPS 
camera on board the Spitzer Space Telescope. We have assembled a sample of 42 clusters from the 
Red-Sequence Cluster Survey-1 over the redshift range 0.3 < z < 1.0 and spanning an approximate 
range in mass of 10 14 ~ 15 M Q . We statistically measure the number of IR-luminous galaxies in clusters 
above a fixed inferred IR luminosity of 2x 10 11 M , assuming a star forming galaxy template, per unit 
cluster mass and find it increases to higher redshift. Fitting a simple power-law we measure evolution 
of (1 + z ) 51±1 - 9 over the range 0.3 < z < 1.0. These results are tied to the adoption of a single 
star forming galaxy template; the presence of AGN, and an evolution in their relative contribution 
to the mid-IR galaxy emission, will alter the overall number counts per cluster and their rate of 
evolution. Under the star formation assumption we infer the approximate total SFR per unit cluster 
mass (SSFR/M c i ustcr ). The evolution is similar, with SSFR/M c i uster - (1 + z) 54±19 . We show that 
this can be accounted for by the evolution of the IR-bright field population over the same redshift 
range; that is, the evolution can be attributed entirely to the change in the in-falling field galaxy 
population. We show that the £SFR/M c i ustCT (binned over all redshift) decreases with increasing 
cluster mass with a slope (£SFR/M c iuster ~ M^ste, A ) consistent with the dependence of the stellar- 
to-total mass per unit cluster mass seen locally. The inferred star formation seen here could produce 
~5-10% of the total stellar mass in massive clusters at z = 0, but we cannot constrain the descendant 
population, nor how rapidly the star-formation must shut-down once the galaxies have entered the 
cluster environment. Finally, we show a clear decrease in the number of IR-bright galaxies per unit 
optical galaxy in the cluster cores, confirming star formation continues to avoid the highest density 
regions of the universe at z ~ 0.75 (the average redshift of the high-redshift clusters). While several 
previous studies appear to show enhanced star formation in high-redshift clusters relative to the field 
we note that these papers have not accounted for the overall increase in galaxy or dark matter density 
at the location of clusters. Once this is done, clusters at z ~ 0.75 have the same or less star formation 
per unit mass or galaxy as the field. 

Subject headings: Galaxies: clusters: general - Galaxies: evolution - Galaxies: starburst - Infrared: 
galaxies 



1. INTRODUCTION 

It is now clear that the environment is a primary factor 
in galaxy evolution, either through direct influences, or 
because galaxy density is a tracer of the important under- 
lying drivers of evolution, such as galaxy mass or forma- 
tion time (e.g., iGomez et all 12003b iBaldry et alJl2006f l. 
Moreover, the importance of the environment may be 
a strong function of cosmic epoch; even if environmen- 
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tal d ependencies are constant with redshift (jPeng et al.l 
I2010T ). galaxies themselves are located in progressively 
more dense regions with time. On the other hand, prop- 
erties of galaxies such as their mass are also predictors 
of their evolution and trends with density are in part 
driven by underlying ma ss-density biases. (jPeng et al.l 
l2010t iMuzzin et all 120121 ) . Understanding the complex 
interdependencies of mass and environment over the his- 
tory of the universe, and their effect on galaxy formation, 
is an immense observational endeavor, requiring substan- 
tial dynamic range in galaxy and halo mass, environment 
or density, and time. 

At z = 0, s tar formation is suppres sed in high-density 
regions (e.g.. iKauffmann et al.l 120041 ) . but recent stud- 
ies have shown that by z ~ 1 star formation has be- 
gun to migrate from low to higher- den sity environments 
(jElbaz et all [20071 ICooper et al.ll200cl . Still, the high- 
est density regions at any given epoch - galaxy clusters 
- are difficult to probe precisely because these regions 
are so rare. While wide-field surveys probe several or- 
ders of magnitude in galaxy density, most do not contain 
large numbers of galaxy clusters, and for these densi- 
ties targeted investigations are more efficient. However, 
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the assessment of the star formation rates of galaxies in 
high-redshift clusters has been hindered in the past by 
the inhomogeneity and sparseness of cluster samples, and 
the different approaches used to quantify the amount of 
star formation (SF) occurring in the cluster environment. 
Conclusions drawn from the comparison of small num- 
bers of clusters selected through different biases are prob- 
lematic. A number of lines of evidence suggest that the 
star formation efficiency of a cluster is strongly correlated 
with its total halo mass (though not necessarily driven 
by it) (e.g.. iPoggianti et al.l 2006). The dynamical state 
of a cluster at the epoch of observation is likely also 
very important, as a major merger event could have a 
profound though temporary effect on the member galax- 
ies. Comparisons between different studies are further 
complicated by star formation rates that are estimated 
through different observational diagnostics such as opti- 
cal emission lines ([Oil], Ha) or infrared emission. These 
different techniques are biased to different galaxy popu- 
lations, or may yield different values for the same galaxy 
either through errors in the calibrations or by probing 
different star-forming regions. Finally, studies have em- 
ployed a variety of different cluster member identifica- 
tion methodologies: spectroscopic confirmation, photo- 
metric redshifts or color selection, each with different 
observational limits and completeness functions. These 
different approaches have all shown some level of evolu- 
tion in the SFR of clusters to higher redshift but have 
varied significant ly in their a ssessmen t of its m agnitude: 
f.gg oc (1+z) 2 - 7 (IKodama et alll200ilGeach et al.| [20F 



Saintonge etaI1l20q8HBai et al.ll2009TlHaines etaLlbooa 
Popesso et al.ll2012t) 

Here we present a study of the average IR properties 
of a large sample of galaxy clusters (42) drawn in a sys- 
tematic way from the first Red-Sequence Cluster Survey 
(jGladders fc Yeell2005[ ) over the redshift range 0.3 < z < 
1.0. The intent of this work is to undertake a very sim- 
ple analysis of the global evolution of dust enshrouded 
activity in cluster environments with time. Our adopted 
method offers a number of advantages, but is also limited 
in scope and for clarity we briefly summarize it here. We 
perform a simple statistical measurement of the back- 
ground subtracted IR counts along the line of sight to 
galaxy clusters. This provides a clean measurement of 
the dust-enshrouded activity in clusters that does not 
require any selection or assumptions beyond a simple 
IR-luminosity cut. We do not, for example, require an 
optical counterpart for the IR sources, which would bias 
the sample toward less dust-enshrouded objects; nor do 
we directly isolate cluster members, which then requires 
a member completeness correction and is again biased. 
Thus, we will be sensitive to all activity above our detec- 
tion limit. Throughout we make the simplification that 
the IR luminosity is entirely produced by star-formation, 
with no AGN conta mination. This assum ption is likely 
an oversimplification (jTomczak et al.ll2011[ ) and does not 
affect the primary conclusions of the paper which depend 
only on the IR luminosity, not the source of activity. We 
discuss AGN throughout the paper, when they are rele- 
vant to the conclusions. 

On the other hand, this statistical method limits us 
to the cluster averaged activity and cannot tell us any- 
thing about the individual galaxies. Therefore we cannot 
control for galaxy specific properties such as stellar mass 



and this limits the implication of the results. More de- 
tailed analysis of the IR galaxy population in this cluster 
sample will be undertaken in future papers, nevertheless, 
important conclusions may be drawn from this simple 
study alone. 

The paper is laid-out as follows. In §2 we outline the 
RCS cluster sample selection and describe the Spitzer 
observations, data reduction and source extraction. In 
§3 we describe the number count analysis methodology 
and in §4 we present the results. §4.1 presents the evo- 
lution of the statistical excess of IR galaxies seen in the 
RCS cluster fields and describes sources of systematic er- 
ror in this analysis; §4.2 converts these measurements to 
the evolution of the integrated star formation per unit 
cluster mass; §4.3 compares this evolution to that of the 
IR-bright field population; in §4.4 we look at the depen- 
dence of these results on cluster richness or mass; and 
§4.5 presents the radial distribution of the IR-luminous 
population in clusters. In §5 we discuss the implications 
of these results and in §6 we summarize the primary con- 
clusions of the paper. We use H = 70 km/s/Mpc, Om 
= 0.3 and il\ = 0.7 throughout. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. The Red-Sequence Cluster Survey Sample 

The cluster sample was drawn from the first Red- 
Sequence Cluster Survey (RCS-1). The RCS-1 is a 90 
square degree optical imaging survey conducted in two 
filters (R c and z' ) with the CTIO and CFHT telescopes 
(jGladders fc Yed I2000D . It was designed to optimize 
the detection of galaxy clusters in the redshift range 
0.3 < z < 1.2 through the detection of the red sequence 
of early-type galaxies within the cluster core. The cluster 
finding method has proven to be extremely robust and 
suffers from minimal projection effects; X-ray imaging, 
spectroscopic verificati on and simulations consistently 
show this to be < 10% ([Hicks et alj|2008t iGilbank et at] 
2007). Moreover, the localization of the red-sequence in 
color space constrains the cluster redshift to ~5%. The 
two-band photometry further provid es an estimate o f the 
richness through the parameter N re d ()Lu et al.H2009T ) , the 
number of red-sequence galaxies brighter than M*+ 2 
within 0.5Mpc. Throughout the paper we use the RCS-1 
estimated N rc( j which uses the z' magnitude and (i? c - 
z 1 ) color to define red sequence galaxies. We note that 
this is a version of the cluster catalog updated from pre- 
vious publications by our group, and will be presented 
in Barrientos et al. (in preparation). The richness mea- 
surements will be discussed in detail in Ellingson et al. 
(in preparation) and Gilbank et al. (in preparation^. 

From the RCS-1 parent sample, we selected 42 clus- 
ters to uniformly fill the redshift range 0.3 < z < 1.0 and 
to span a richness range which corresponds to roughly 
an order of magnitude in mass (N rc d =10-60, and us- 
ing the N ro d to mass conversion outlined below, M ~ 
10 14 ~ 15 Mq, with ~ 30% uncertainties). In so doing, we 
have assembled a sample which is no longer representa- 
tive of the RCS selection distribution in mass or redshift, 

1 We note that this richness measurement is essentially the same 
as the Bgc parameter we have used in previous work but uses a 
somewhat larger counting radius and does not extrapolate the lu- 
minosity function beyond the observational limits. Our previous 
Bgc measurements may be approximately converted to N re( j by 
dividing the former by ~30 



but rather comprises a representative sample with the 
goal of isolating the effects of redshift and cluster prop- 
erties on cluster galaxy evolution. The size of the sample 
was motivated by the desire to remove uncertainties due 
to cluster-to-cluster variations, by averaging over bins of 
redshift and/or richness with a minimum of ten clusters 
per bin. Figure Q] shows the richness and redshift values 
for the cluster sample; no selection bias in richness with 
redshift is evident. Table [1] summarizes relevant cluster 
properties. 

The RCS-1 clusters are the focus of a number of spec- 
troscopic campaigns by the RCS consortium for the pur- 
pose of redshift confirmation, population and cluster dy- 
namic studies, and gravitational lensing analyses, which 
have provided extensive but inhomogeneous coverage of 
the RCS-M1PS fields. In Figure [I] and Table [I] we indi- 
cate the 25/42 clusters with spectroscopic redshift confir- 
mation. The larger spectroscopic sample have also cal- 
ibrated the relation between the N ro( j richness and the 
velocity dispersion and this is presented in Ellingson et 
al. (in preparation). 

In this paper we work within a radial limit of r2oo of 
each cluster and also normalize by the cluster M200 mass. 
The effects the uncertainties on these properties have on 
our results are discussed in detail in the Appendix. We 
calculate the r2oo radii by first estimating the velocity 
dispersion from the optical richness as outlined above. 
The r2oo radii were determined through the following re- 
lation (|Carlberget alJll996t l: 
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2.2. Spitzer Imaging: MIPS 24fJ-rn and IRAC 
3.6/4.5fim 

Our primary data set, the Spitzer 24/im imaging, 
was obtained through open time program 30940. The 
MIPS observations were designed to reach an approxi- 
mate depth of Lir ^lO 11 L© out to a radius of r2oo for 
each cluster and therefore per pixel integration times and 
image sizes vary depending on cluster richness, redshift 
and the thermal background. A subset of the clusters ob- 
tained longer integrations to facilitate a deeper study on 
a smaller number of clusters, but this additional depth 
is not relevant for this work. The integration time and 
map area for each cluster are listed in Table [1] aver- 
age exposure times per pixel range from 200s to 4750s 
and image areas from 50 arcmin 2 to 300 arcmin 2 , total- 
ing ~1.5 deg 2 . The MIPS images were reduced using 
a combination of the Spitzer Science Center's MOPEX 
software and our own IDL routines developed to further 
optimize background subtraction in each field. 

In §4.4 and the Appendix we will compare the MIPS 
measurements to a comparison IRAC-selected (rest- 
frame NIR) population. For this we use deep IRAC 4- 
channel imaging which was obtained through open time 
program ID 20754 for all z > 0.5 clusters within the RCS- 
MIPS sample, with the exception of RCS022056 and 
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Fig . 1. — The optical richness, parametrized by N rec j IILu et al.l 
2009) versus the redshift of each cluster in the RCS-MIPS sam- 
ple presented here. Solid points denote clusters with spectroscopic 
redshift confirmation and open points correspond to redshifts de- 
termined through a fit to the color location of the red sequence. 
The corresponding M200 was computed following the N re[ j-cr rela- 
tion of Ellingson et al. (in preparation). The uncertainties can be 
applied to both axes, but in the case of M200 do not include any 
intrinsic scatter in the richness-mass relation. 



RCS022158. These two clusters are located within the 
public SWIRE fields and we do not use these two clusters 
in the IRAC-dependent analysis to avoid introducing sys- 
tematic errors due to the different observational mode of 
the SWIRE data. The goal of the program was to probe 
the stellar mass of the cluster galaxies and therefore the 
IRAC sample was limited to z > 0.5; at lower redshifts 
NIR imaging samples the appropriate rest-frame wave- 
length and these clusters are part of an on-going CFHT 
imaging program. The IRAC and CFHT observations 
were designed to reach an approximate depth of M*+1.5 
at each redshift. 

Clusters with redshifts below z ~ 0.78 were centred 
on the 3.6^m array and those above this redshift were 
placed on the 4.5^m array; this ensured sampling at the 
observed wavelength which most closely corresponded 
to rest-fra me K. Image processi ng; was performed using 
IRACproc (jSchuster et al.l [20 061 a software suite devel- 
oped to wrap the existing MOPEX pipeline in IRAC 
mode and add IRAC specific improvements to the outlier 
rejection. 

2.3. Source Detection and Photometry Catalog 

The 2AfjLTn source extraction algorithm combined the 
source-finding capabili ties of the PPP detection and pho- 
tometry program (Ycc 1991) with the aperture photom- 
etry and PSF-fitting capabilities of DAOPHOT (|Stetsonl 
|1987[) . Source detection and apertur e photome try on the 
IRAC images was done using PPP (Yce 1991) alone. 

Completeness limits at 24//m were estimated through 
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TABLE 1 

Information on cluster sample and Spitzer-MIPS observations 



Cluster ID 


rcdshift a 


N red 


?"200 


integration 


map size 








(arcminutes) 


time (s) 


(arcminutes 2 ) 


RCS212134-6335.8 


0.217 


27.5 ± 7.3 


7.6 ± 1.8 


200 


225 


RCS035139-0956.4 


0.304 


23.3 ± 6.6 


5.2 ± 1.3 


200 


225 


RCS022516+0011.5 


0.357 


20.6 ± 6.6 


4.3 ± 1.2 


200 


225 


RCS132655+3021.1 


0.35 


23.6 ± 6.8 


4.6 ± 1.2 


200 


300 


RCS144726+0828.3 


0.376 


50.3 ± 9.4 


6.5 ± 1.1 


200 


300 


RCS092821+3646.5 


0.393 


27.0 ± 7.8 


4.9 ± 1.3 


200 


300 


RCS022359+0126.1 


0.394 


19.0 ± 6.2 


3.7 ± 1.1 


200 


225 


RCS145226+0834.6 


0.395 


14.8 ± 5.4 


3.3 ± 1.1 


200 


225 


RCS051834-4325.1 


0.396 


21.8 ± 7.0 


4.1 ± 1.2 


200 


300 


RCS022403-0227.7 


0.408 


23.6 ± 6.8 


4.1 ± 1.1 


200 


225 


RCS231526-0046.7 


0.40 


40.1 ± 8.5 


5.4 ± 1.0 


200 


300 


RCS044207-2815.0 


0.45 


29.7 ± 7.4 


4.4 ± 1.0 


200 


100 


RCS215223-0503.8 


0.480 


51.8 ± 9.6 


5.4 ± 0.9 


200 


150 


RCS051855-4315.0 


0.508 


21.8 ± 6.5 


3.3 ± 0.9 


200 


100 


RCS110733-0520.6 


0.511 


17.5 ± 6.0 


2.9 ± 0.9 


600 


100 


RCS234717-3634.4 


0.55 


50.4 ± 2.6 


9.6 ± 0.4 


200 


150 


RCS110104-0351.3 


0.571 


19.5 ± 6.3 


2.7 ± 0.8 


600 


150 


RCS144654+0827.0 


0.628 


17.5 ± 5.8 


2.4 ± 0.7 


200 


100 


RCS144557+0840.3 


0.629 


29.0 ± 7.2 


3.2 ± 0.7 


200 


150 


RCS110439-0445.0 


0.637 


24.4 ± 6.9 


2.9 ± 0.7 


600 


100 


RCS215248-0609.4 


0.649 


22.6 ± 6.4 


2.7 ± 0.7 


3300 


100 


RCS211852-6334.6 


0.658 


21.4 ± 6.9 


2.6 ± 0.7 


200 


100 


RCS110246-0426.9 


0.70 


15.3 ± 5.7 


2.1 ± 0.7 


600 


200 


RCS212238-6146.1 


0.70 


23.4 ± 7.1 


2.6 ± 0.7 


200 


100 


RCS112225+2422.9 


0.70 


23.0 ± 6.5 


2.5 ± 0.6 


2000 


100 


RCS141910+5326.1 


0.710 


31.8 ± 7.6 


2.5 ± 0.4 


2000 


100 


RCS234220-3534.3 


0.70 


34.0 ± 8.0 


3.0 ± 0.6 


200 


100 


RCS044126-2813.2 


0.734 


21.3 ± 6.6 


2.4 ± 0.7 


100 


100 


RCS132939+2853.3 


0.75 


13.6 ± 2.0 


5.2 ± 0.7 


200 


50 


RCS022433-0002.3 


0.773 


24.6 ± 6.8 


2.0 ± 0.5 


3300 


100 


RCS110411-0337.5 


0.80 


22.7 ± 6.7 


2.4 ± 0.6 


600 


100 


RCS051940-4402.1 


0.827 


20.5 ± 6.8 


2.1 ± 0.6 


200 


50 


RCS110118-0328.6 


0.80 


22.0 ± 6.7 


2.2 ± 0.6 


600 


50 


RCS110206-0414.5 


0.90 


15.5 ± 5.7 


1.8 ± 0.6 


600 


50 


RCS110615-0330.8 


0.90 


20.0 ± 6.3 


2.0 ± 0.6 


600 


100 


RCS162009+2929.4 


0.869 


23.1 ± 6.7 


2.2 ± 0.6 


2000 


100 


RCS231953+0038.0 


0.907 


39.4 ± 8.4 


2.8 ± 0.5 


3200 


100 


RCS132631+2903.3 


0.919 


28.1 ± 7.2 


2.3 ± 0.5 


200 


100 


RCS022158-0340.1 


0.90 


25.5 ± 7.1 


2.1 ± 0.5 


4750 


50 


RCS051908-4323.3 


1.0 


21.1 ± 6.8 


1.8 ± 0.5 


200 


100 


RCS022056-0333.3 


1.0 


29.2 ± 7.5 


2.2 ± 0.5 


4750 


50 


RCS043938-2904.8 


0.956 


26.5 ± 7.1 


2.2 ± 0.5 


2000 


100 



a Redshifts in bold are spectroscopically confirmed; non-bold values denote redshifts deter- 
mined through a fit to the location of the cluster red-sequence, good to ~5%. 



the insertion and recovery of fake sources into each clus- 
ter image, using the same detection and photometry 
method as used for the un-altered images. An input 
source was deemed recovered if an object was found 
within 2" of the input position. No secondary check 
against the true source catalog was made, and no in- 
put/output flux ratio criterion was employed. While this 
may lead to a slight underestimate of the actual comple- 
tion limits it requires the fewest assumptions and subjec- 
tive definitions of a recovered object. Figure [2] shows the 
80% completeness limit for each cluster field, compared 
with the expected 24/im flux of an infrared-luminous 
galaxy with Ltr ~ 2xlO n L^ over the same redshift 
range (jCharv fc ElbazJl200H iDale fc Heloull2002h . 

Fields with two MIPS AORs taken at different times 
reveal a small number of asteroids. Theoretical pre- 
dictions of the Tedesco' s Statistical Asteroid Model 
(jTedesco &: Zap pala 20 051) and the results of other ex- 
tragalactic programs (jPapovich et al.ll2004T ) indicate that 
the number of asteroids per field as well as the spread in 
number density between the ecliptic and zenith does not 



significantly alter the galaxy counts. We therefore do not 
attempt to identify or remove the asteroid contribution 
in the cluster fields with multiple AORs. 

2.4. Background Comparison Field: SCOSMOS 

The analysis presented here requires measurements of 
the background field galaxy population (i.e. along the 
line of sight). We use the p ublicly available SCOSMOS 
G03 24/zm data ([Sanders et al. 2007) which is of com- 
parable depth as our deepest cluster images and cov- 
ers 2 deg 2 (Figure [3]). To avoid introducing system- 
atic errors into our analysis we worked directly with 
the original SCOSMOS unreduced data (rather than the 
publicly available source catalogs or images) and re- 
duced/mosaicked them in an identical manner as the 
cluster fields. 

The reduced SCOSMOS comparison field was then 
passed through the same detection and photometry 
pipeline as the clusters. In Figure [3] we show the result- 
ing differential number counts (Euclidian-normalized) of 
the cluster fields and our SCOSMOS catalogs and com- 
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Fig. 2. — An illustration of the MIPS image depth for each cluster 
field compared to the expected 24fira flux of a template galaxy with 
Ljr = 2xl0 11 Lp 1 ; the solid line correspond s to t he co mbined mod- 
els of Chary & Elbaz (2001) and Dale & Hclou (2002) and the two 
dotted lines show the range in these two model libraries. The solid 
points show the ~80% completeness level for each separate clus- 
ter field as determined by the Monte Carlo simulations described 
in the text. The dashed line denotes the Lj R parametrized depth 
used in §4.2: we set the z = 1 depth to Ljr = 2x10 Lq and alter 
the z < 1 depths to follow the evolution of Ljn of the field. In 
§4.1 we count all galaxies above the solid line and in §4.2 we count 
above the dashed line. 



pare with the published counts of lPapovich et all ((2004). 
Good agreement is seen between our SCOSMOS analysis 
and Papovich et al., while the clusters show a slight over- 
all excess which will be discussed throughout the paper. 
We have also averaged the cluster counts over four expo- 
sure time bins to again illustrate the difference in depth 
within the cluster sample. The flux at which the cluster 
counts diverge from the field counts due to completeness 
(i.e. where the cluster count curve falls below the field) 
for each bin agrees well with the MC results described 
above. 

3. ANALYSIS: GALAXY COUNTS AND BACKGROUND 
SUBTRACTION 

To estimate the observed 24/im flux limit at each 
cluster redshift, corresponding to a constant luminos- 
ity limit, we e mploy the models and pr e script ions of 
iCharvfc Elba3 ((2001 and iDale fc Heloul pOOl . Im- 
plicit in this is the assumption that the shape of the 
lR SED of galaxies does not evolve appreciably over the 
redshift range of interest. We work to the limit of Ltr 
~2xlO n L Q : this is set by the completeness limits of 
the cluster MIPS images (Figure [5]) and indeed the ex- 
act luminosity limit is less importa nt than its c onsistency 
across all cluster fields. Following iBelll (|2003l ) this limit 
corresponds to a star formation rate of ~30 M Q yr _1 ; 
although rare in the local universe this is a more repre- 
sentative level of field galaxy star formation by z > 0.5 



Fig. 3. — The Euclidian- normalized differential number c ounts 
of MI PS sources in the cluster fields compared to ( Papovich et al. 
[2001) (blue diamonds) and our analysis of the SCOSMOS G03 
data (solid orange circles). The uncertainties on the Papovich 
counts do not include field-to-field variance. The open circles cor- 
respond to the cluster fields divided into four integration time bins: 
200s, 600s, 2000s, > 3200s. The size of the circle increases with in- 
tegration time. The two solid lines connect the counts for the 200s 
and 600s images, to draw the eye to the region where they begin 
to diverge from the SCOSMOS counts due to incompleteness, at 
~100^iJy. The longer integration time images are consistent with 
COSMOS to lower levels. Also shown (dotted lines) are the in- 
ferred luminosities of a galaxy with Lir = 2xlO 11 L0 (LIRG) for 
z = 0.6, 0.8 and 1.0 that we employ in this work (2 = 0.4 lies off 
the plot). We count galaxies to the right of the dotted lines for 
each redshift bin, and thus this illustrates the sufficient depth of 
all of the images, given the adopted limit, and the adequate depth 
of the comparison SCOSMOS field. 



(|Le Floc'h et alJl2005[ ) . In this analysis we have assumed 
the contamination of the MIR emission from AGN is neg- 
ligible. In doing so we may be introducing systematic un- 
certainties into the measurements and this is discussed 
in more detail in the Appendix. 

Background counts are determined in a similar man- 
ner. For each cluster we measure the average number 
of galaxies above its calculated flux limit and within the 
corresponding radius in 500 randomly placed apertures 
in the COSMOS field. The standard deviation of these 
counts provide us with an estimate of the field-to-field 
variance in the background, and this is what we take to 
be the dominant uncertainty on each cluster measure- 
ment. We do not attempt to incorporate the error in 
^200 stemming from the richness uncertainties, the in- 
trinsic scatter in the mass-richness relation, any random 
error in the redshift, nor any systematic error in the flux- 
limit estimate. Random uncertainties will be minimized 
by our stacking technique and systematics only matter if 
they are redshift dependent. We discuss the sources of 
error in detail in Appendix A. 

4. RESULTS 
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4.1. The number counts of IR-Luminous galaxies in 
cluster fields 

In Figure |4] we show the total number of galaxies with 
Lir > 2x 10 11 Lq in excess of the average background 
counts, within r 2 oo of each cluster center. We choose r2oo 
as the largest radius for which all clusters have uniform 
MIPS coverage. To remove additional scatter due to the 
random variation in the average mass from bin to bin 
we further normalize each integrated count by the M200 
mass of the cluster (inferred from the richness as outlined 
in §2.1). 

The top panel shows the measurements of each indi- 
vidual galaxy cluster and the scatter from cluster to clus- 
ter is substantial; given the size of the error bars this is 
dominated by the variance in the background, though 
real differences between the clusters must also be impor- 
tant. The increase in scatter with redshift is expected 
as a fixed luminosity limit corresponds to a deeper flux 
limit at higher redshift and the angular size correspond- 
ing to a given r2oo decreases, therefore the shot noise in 
the background increases. The blue region denotes the 
68% range in measurements conducted on random areas 
of the COSMOS image. In spite of the limitations inher- 
ent in statistical background subtraction, a clear increase 
in the average number of IR luminous galaxies per unit 
cluster mass is seen with redshift (bottom panel). To 
quantify the evolution we fit a simple power-law, N (> 
S) = N (1+z)", and find n = 5.1±1.9, where the un- 
certainties correspond to the la significance region in a 
numerical chi-squared fit. We also show the resulting fit 
if the cluster with the highest value within each bin is re- 
moved (i.e. 10% of the sample), as the lower dashed line. 
In this case the evolution goes as n =4.8±2.0, indicat- 
ing that the trend is not due to a single outlier in each 
bin. Finally, we note that the chosen power-law func- 
tion is arbitrary, and though it allows a quantification of 
the average evolution smoothed over redshift that can be 
compared to other work, it should not be interpreted as 
the evolution of a single cluster. Indeed, as mentioned 
below, the true evolution may be much more stochastic 
as, for example, groups are accreted onto clusters. To 
further illustrate this we simply show the data in two 
redshift bins (open points) where the difference between 
high and low redshift (divided at z = 0.7) is 2.5a. 

The absolute counts of IR galaxies and the rate of evo- 
lution, is dependent on the adopted SED. As outlined 
above, we have used the now-standard libraries of Chary 
& Elbaz and Dale & Helou, but these do not include 
AGN contributions to the MIR emission. If there is sig- 
nificant AGN contamination then we are in fact measur- 
ing the number counts to a different intrinsic flux level 
(lower Lir for a given observed 24/^m flux) , and that flux 
level may not be constant with redshift, in particular if 
the rate of AGN contamination also evolves. This is ex- 
plored in more detail in the Appendix, where we conclude 
that we do not see evidence for variable contributions 
of AGN in our sample - though the statistics are poor. 
Still, several studies have indicated that AGN are not a 
major contaminant over t he redshift and flux range ex- 
plored h ere. For example. iBrand et al.l ([2006D and more 
recently iKirkpatrick et al.l (|2013t ) conclude that at the 
24^m flux levels of interest here (200-600/iJy), signifi- 
cant AGN are present in only 10% of galaxies. Brand 



et al. also estimate the amount of contamination AGN 
contribute to the total MIR flux to be ~10% for similar 
depths. Finally, we note that because our technique re- 
lies on background subtraction of the statistical field, sys- 
tematic differences between the two (as seen here) would 
require corresponding differences in the IR SEDs of field 
and cluster galaxies - and would therefore remain an in- 
teresting result - though requiring a different interpreta- 
tion. Understanding the importance of AGN is extremely 
important and will be addressed in later work using more 
appropriate data than presented here. 

For the remainder of the paper, we assume that the 
MIR emission is produced entirely by star formation. 
In doing so, these measurements can further provide an 
estimate of the evolution of the more physically mean- 
ingful quantity, the total SFR per unit cluster mass, 
SSFR/M c i ustcr and this is shown in Figure [5]). To 
estimate this we use the number counts of Figure U] 
to set the normalizati on of the IR luminosity function 
(|Le Floc'h et al.l I2005D at each redshift. We then esti- 
mate the total Lir from galaxies above 2x 10 11 L Q by in- 
tegrating the luminosity function. We then convert this 
to a ESFR following Bell (2003). We again fit a simple 
power-law and find roughly the same exponent as before: 
n = 5.4 ±1.5. 

4.2. Comparison with the Field 

We see a clear increase in the average number of IR- 
luminous galaxies in clusters with redshift, but this evo- 
lution must be compared with the behavior of field galax- 
ies at similar redshifts. Since clusters are continually 
accreting from the field any change in their galaxy pop- 
ulations with time may simply reflect the corresponding 
evolution in field galaxies feeding the infall, rather than 
cluster specific evolution. 

To look for differences between field and cluster pop- 
ulation evolution we use the s a me m ethod of the lower 
redshift study of IHaines et al.l (|2009Q . Instead of mea- 
suring the number of galaxies within clusters to a fixed 
luminosity at all redshifts, we vary the intrinsic depth 
to reflect the evolutio n in the IR luminosity function of 
ILe Floc'h et~aH (|2005l ): L* R ~ (1 + z) 32 . We normalize 
to the previous depth of Lir — 2x 10 11 L Q at z — 1 
and show the corresponding flux depth with redshift in 
Figure This has the added advantage that it allows 
us to reach the intrinsically deeper luminosity limits in 
the lower redshift images. The results of this analysis 
are shown in Figure [5] (upper black points). An excess of 
IR bright galaxies is seen in all redshift bins, reflecting 
the over-density of galaxies at the location of the clus- 
ters, but there is now no evidence for an increase in the 
number with redshift. A formal fit of N = N (l+z) n 
to the data yields n — 0.5±1.4, indicating the cluster IR 
population follows the field evolution. 

For further comparison, we estimate the number of IR 
galaxies per unit 10 14 Mm in the field by int egrating the 
luminosity function of lLe Floc'h et all (J2005I ) at each red- 
shift, and normalizing by the average matter density of 
the universe. Note that this is different than the back- 
ground subtraction which removes galaxies along the line 
of sight, and the correction form the field IR luminos- 
ity function, which removes overall expected evolution of 
galaxies. This simply scales the star formation per unit 
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Fig. 4. — The mean number of objects within T200 in excess of the background counts with Lir > 2x10 Lq. The upper panel shows 
the measurements for each individual cluster. Error bars represent the rms variation of the field counts and the blue band encompasses the 
68% region of the same analysis done on random locations in the COSMOS field. The scatter in the mass-richness relation (see Appendix 
A) is not included here. The bottom panel shows the weighted average of the clusters in four redshift bins (solid points) and within two 
redshift bins (small open points), with the uncertainties corresponding to the standard error in the mean. For the latter case the low and 
high redshift points are different at a l evel o f 2.5<r. The top dotted line corresponds to the evolu tion of the fracti o n of s tar-forming galaxies 
in clusters measured by Haines ct al. (2009) at z < 0.3 using LoCuSS clusters and the data of Saintongc ct al. (2008), normalized to our 
measured value at z = 0.4. The middle black line shows the best fit power-law to the four RCS bins, of the form N(> S) ~ (l+z) n : n = 
5.1 , and this is also shown in the upper-panel (dotted line). The bottom dashed line shows the trend when the highest cluster in each 
bin is removed from the analysis. The two stars at z ~ 1 denote the results of the extinction-corrected optical study of GCLASS clusters 
(Muzzin ct al. 2012): the lower star corresponds to the confirmed number of star forming galaxies above 30 M0yr — 1 (using the [Oil] line 
measurements) and the upper star shows the value corrected for spectroscopic completeness (see discussion in text for more details). The 
dotted line in the upper panel is identical to the solid line (best-fit function) in the lower panel. 
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Fig. 5. — The total star formation rate per unit cluster mass, 
£(SFR)/M, shown averaged over all clusters in four redshift bins. 
The solid line shows a fit to these data of the form SSFR/M <x 
(1 + z) n with the best fit value n = 5,4± 1.5. The dot-dashed 
line shows the results of the similar study by Podcsso ct al. (2012) 
using 9 clusters over the same approximate redshift range, and 
the dashed line is the same scaled to ou r IR depth and aver age 
cluster mass. The two stars correspond to[Muzzin ct al. (2 0121 '). as 
in Figure |4] 
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Fig. 6. — Similar to Figure U but with the evolution of the IR 
luminosity function derived from field galaxies accounted for in the 
IR count depth. We measure the number of IR luminous galaxies 
to the same depth relative to an evolving L^. The black solid 
line shows the best fit power law: N~(l+z) n to the data yields 
n = 0.5±1.4. The blue points show the cluster values after the 
mean number of IR-luminous galaxies per halo mass in the field 
has beens subtracted; see 54.2 for details. 



mass in the field to the mass of the cluster. The difference 
between this field estimate and the cluster measurements 
is shown in blue and does not indicate an overall excess 
or deficit in the number of IR bright galaxies in clus- 
ters, once normalized for the halo mass. Note, however, 
that while we have normalized the cluster counts to the 
mass within r2oo we may actually be sampling a larger 
mass along the line of sight, by up to 30%, if the clus ters 
are not highly concentrated ( Lokas fc Ma mon 200T|). If 
so, we have then under-subtracted the field and the blue 
points should be lower, by roughly the same percentage, 
possibly indicating a small deficit of IR galaxies within 
clusters. While severa l studies (e.g. iKocevski et al.ll20lH 
iMarcillac et al.l I2007T) claim to measure enhanced star 
formation in z ~1 clusters, relative to the field, they 
have neglected the crucial step of accounting for the over- 
all increase in mass (stellar or dark matter) in clusters. 
Here we show that the overall star formation per unit 
halo mass, or per unit galaxy density, is either consistent 
with or lower than the field value; that is, any apparent 
increase in star formation within this sample of clusters 
at 0.5 < z < 1.0 is accounted for simply by the overall 
increase in galaxies at that location. 

4.3. Correlation with Cluster Richness/Mass 

In Figure [7] we show the integrated star formation rate 
per unit cluster mass, in four mass bins. As before, this 
includes all objects with Lir > 2xf0 11 L Q and within 
?*200> with the background subtracted. We have further 
normalized all counts to the counts of z = 0.4 to remove 
the redshift evolution, as we do not have enough clus- 
ters to bin in redshift and richness, and still maintain 
decent statistics. This last step is not strictly necessary 
provided all mass bins represent the same redshift dis- 
tribution, but doing so should reduce the scatter due to 
small numbers. We use dynamic bin sizes to ensure equal 
numbers of clusters in each bin; as the mass distribution 
is not uniform over the range, this results in a larger bin 
width for the highest masses. We fit a simple power-law 
and find £(SFR)/M cl ex M^* " 4 . 

As discussed in the Appendix, systematics in the mass 
normalization could bias the total SFR per unit cluster 
mass. In Figure [7] however, systematic errors are not 
required to produce a trend; a false decrease in the inte- 
grated SFR per unit mass with cluster mass could be in- 
troduced through statistical uncertainties or through the 
intrinsic scatter (~30-40%) in the ma ss-richness relation 
(|Rozoet all 120 lit iRvkoff et al1l2012f ). Consider a sam- 
ple of clusters of equal true mass but with a 30-40% ran- 
dom error on their richness measurements. This would 
spread the sample non-uniformly over a larger mass range 
(though not as large as explored here) and when binned 
in mass this results in a systematic over estimate of the 
mass in the high-mass bins, and a systematic underes- 
timate of the mass in the low-mass bins. As discussed 
before this will result in an increase of S5Fi?/Af c i ustcr at 
lower masses and a decrease at higher masses, as is seen. 
Still, once again, the measured slope is steeper than what 
can be accounted for by reasonable uncertainties alone. 
For a flat radial distribution of SF galaxies, the expected 
slope relation is S SFR/M c i ust0 r oc M -03 . Reproducing 
the observed relation requires an unreasonably steep ra- 
dial distribution of galaxies (so that the number counts 
are dominated by the inner regions and no new counts 
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Fig. 7. — The integrated star formation rate of clusters per unit 
cluster mass as a function of total cluster mass. As elsewhere, this 
includes all galaxies above Ljr > 2x 10 11 Lq and within r2oo- 
The richness-mass conversion method is described in the text. We 
have removed the effect of redshift evolution by normalizing all 
cluster counts to the z = 0.4 level, following the relation measured 
in §4.1. The dotted line corresponds to the best fit power-law of 
SSFR/M clustcr oc M" 1 - 5 * - 4 . 



are incorporated as the radius increases) and uncertainty 
on the richness measurements of more than a factor of 
two. The decrease in the integrated SFR per unit cluster 
mass with cluster mass therefore appears robust. 

4.4. The Radial Distribution of MIPS Galaxies 

In Figure |8] we show the average number density of 
galaxies above the background in distinct cluster-centric 
annuli. We show the distribution for IR luminous galax- 
ies (yellow) and perform the same analysis on the IRAC 
NIR galaxy population at the same locations in the clus- 
ter fields (blue). Both are normalized to the outermost 
radial bin. The NIR galaxies provide a trace of the under- 
lying galaxy surface density, and roughly correlate with 
stellar mass. Note however, that the actual stellar mass 
density will be a steeper function of radius because the 
average stellar mass per galaxy increases to ward the cen- 
tre of the cluster (e.g., M uzzin et al.ll20i2f) . We do not 
split the sample into redshift bins as the statistics become 
too poor: thus there is the implicit assumption that the 
radial distribution is the same over 0.5< z <1.0, even 
though the total number of IR galaxies evolve. Each 
population is analyzed separately and we do not require 
galaxies to be detected in both IRAC and MIPS. Fig- 
ure E] shows that the MIPS galaxies follow a relatively 
flat radial distribution, with a hint of an increase toward 
smaller radii, reflective of the general increase in galaxy 
density in this region, as illustrated by the steep radial 
distribution of NIR-detected galaxies. The black points 
show the ratio of the number density of these two pop- 
ulations - the number of IR-luminous galaxies per unit 



rest-if-band selected galaxy (to the limits given in the 
caption). 

5. DISCUSSION 

5.1. The evolution of the IR-luminous galaxy population 

in clusters 

Here we present the largest single study which at- 
tempts consistency over the redshift range 0.3 < z < 
1.0 in cluster selection, member identification and SFR 
measurement method. We find the absolute number of 
IR galaxies in clusters, normalized by the parent halo 
mass, or the corresponding approximate SFR/M c i ust0 r 
rises as (1 + z)~ 5 . The primary result of this work is 
that although there is seemingly rapid evolution, such 
behaviour is entirely consistent with the evolution of the 
IR bright field galaxy population over the same redshift 
range. Thus the evolution of the IR galaxy population 
in clusters, and the inferred SFRs within clusters, can 
be attributed solely to the change in the in-falling field 
population. 

Given the variety of approaches and samples outlined 
above, it is difficult to directly compare our results to 
those of other a uthors. The observa t ions of single or a 
few clusters fe.g.lKodama et al.l[200llGeach et alJl200d 
iSaintonge et al.ll2008| ) certainly agree very generally with 
our work, given that they reach different SFR limits, 
measure the total SFR to different radii, and select dif- 
ferent cluster masses; but a meaningful comparison must 
be limited to a few surveys which target larger numbers 
of clusters and use methods which are easily corrected to 
our own. The Haines et al. survey of 30 LoCUSS clusters 
below z = 0.4 (Figure H} is a large sample of systemat- 
ically selected X-ray clusters, of a similar mass range as 
explored here. They undertake a statistical estimate of 
the fraction of star forming galaxies in clusters, relative 
to the optical population, which is analogous though not 
identical to our measurement, but reach 2x deeper in IR 
luminosity. They find a very similar rate of evolution as 
found here, of (1 + z) 5 ' 7 , which implies a continuation of 
the higher redshift evolutionary trend to z ~ 0. 

Our agreement with the recent results of lPopesso et al.l 
(2012), shown in Figure[5j is also heartening. They reach 
2x deeper in SFR (also IR-estimated) and thus their 
SFRs must be corrected upwards (using the LIR lumi- 
nosity function at each redshift). The average cluster 
mass of the Popesso et al. sample is also ^3x higher 
than the RCS-MIPS average mass. As we sho wed in the 
previ ous section (and found by others e.g. (|Bai et al.1 
l2007f l) the ESFR/M c i us ter decreases with cluster mass, 
and therefore the Popesso et al. relation must be scaled 
downward (we use the relation found in §4.4). These two 
effects roughly cancel and preserve the good agreement 
between these two studies. 

Finally, in Figure H we a lso show the SSFR /M of 
the GCLASS z -1 clusters (|Muzzin et al.l l2012h . The 
GCLASS clusters a r e derived from the SpARCS survey 
(jMuzzin et al.ll2009t iWilson et aTll2009t ) and are also se- 
lected through the Red-Sequence technique; they there- 
fore make an ideal higher redshift comparison. The star 
formation is optically estimated through [Oil] lines and 
includes spectroscopically identified cluster galaxies with 
[Oil] (extinction corrected), but we cut it at the same 
SFR limit as the MIPS data: 30 Mayr" 1 . Again the 
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agreement with the RCS number counts excellent (see 
Figure[5]). The Popesso and Muzzin programs differ from 
ours in a fundamental way in that the star forming galax- 
ies must be detected in the optical. In the case of Popesso 
an optical counterpart is required and for the GCLASS 
clusters the SFR is also estimated through the strength 
of the [Oil] line. The agreement in the overall normal- 
ization of these three distinct measurements supports the 
idea that optically biased studies are not missing a sig- 
nificant fraction of the star formation above this limit 
out to z ~ 1. Still, this conclusion assumes that above 
a SFR of 30 M Q yr— 1, optical and IR observations are 
measuring the same population of galaxies, and this may 
not be the case: optical studies will be biased to dust-free 
systems and IR observations to dust-enshrouded objects, 
and each may be missing a signifiant fraction of the other. 

5.2. IR-Luminous Galaxies and Cluster Mass Assembly 



But how important is this in-falling IR population to 
the assembly of the stellar mass in clusters at z ~ 0, and 
can we constrain the timescale over which this activity is 
shut down? To investigate this we take an average mass 
halo in our sample, M200 ~ 4xl0 14 M©, and note that 
it will evolve into a p resent-day halo of mass ^10 15 M Q 
(jWechsler et al.ll2002D . By integrating the SFR of the 
IR-bright population over the redshift range probed here 
(taken from Figure [5]) we estimate that this population 
would add ~4xlO 11 M0 of stars to the cluster. If the 
decline continues with the same form to z = the mass 
added increases to ^8xlO n MQ. Observations of local 
clusters in dicate that the stellar mass is ~1% the total 
halo mass (jAndreon et al]l2010T ). or 10 13 M Q and there- 
fore this simple treatment indicates the in-falling infrared 
bright population below z ~ 1 is responsible for 10% of 
the total stellar mass of such a cluster today. However, 
the results of §4.4 indicate that as a cluster grows in mass 
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its integrated SFR per unit mass will decrease following 
Figure [3 If we account for this additional shutting down 
of star formation, due not to the overall change in the in- 
falling population but to the change in the parent halo 
(again using the halo growth models of iWechsler et all 
(2002)) the assembled stellar mass is reduced by roughly 
a factor of two, or to 5%. 

The descedent population is harder to constrain. The 
only galaxy population which shows strong evidence for 
evolution from z ~ 1 is the faint-end of the re d-sequence 
(IGilbank et alll2008t fPe Lucia etlT]|2008UVuicani et al.l 
I2010D . however the estimated mass added by the IR- 
bright population is comparable to t he total mass of the 
dwarf population in clusters at z ~ (jBildfell et al.ll20T2l) 
and would therefore require stronger observed evolution 
in this population than is currently seen, even with the 
above caveats. However, we note that the uncertainties 
in the mass increase in the bright-end of the red-sequence 
over the redshift range probed here are large enough that 
they can in corporate the mass added by the IR-bright 
population (|Bildfell et al.ll2012D . 

Because of our statistical technique we have no in- 
formation about the average current mass in stars in 
these systems; however the average mass of field LIRGs 
at z ~ 0.5-1 is -0.5xlO n M m (IGenzel et al.1 l20lot 
lElbaz et alJlMll . or sub-M* (|Vulcani et alJl2012ft . To 
restrict growth to the faint-end of the red-sequence (or 
<10 n Mq) the star formation seen here (taken to be 
at least 30 M0yr _1 ) must be quenched within < 2 Gyr. 
This could be accomplished through simple gas depletion 
(iBalogh et all 12003 [Genzel et alj|20lot McCarthy et al.l 
|2008f ) provided the replenishment of the gas reservoir is 
shut down once the galaxy enters the cluster environ- 
ment. Still, longer quenching timescales are easily ac- 
commodated since as noted above this population need 
not be restricted to the faint end of the red-sequence. 
Understanding the fate and more detailed importance of 
the IR-luminous population in clusters will be explored 
in future work, with the aid of spectroscopy, for now we 
simply conclude that our results are in line with other 
evolutionary studies of cluster galaxy populations. 

5.3. The dependence of star formation on cluster mass 

The decrease in ESFR/M c i uster with halo mass seen in 
Figure [7] is not surprising: this beh avior has been seen 
by o thers (IFinn et al.ll200l l200l IBai et all [2007) (but 
c.f. IChung et al.l ([2012D ) and indeed is expected from 
the properties of local (i.e. z =0) galaxy clusters. The 
stellar mass frac tion of clusters decre ases with increas- 
ing cluster mass ()Andreon et alll2010| ) as does the metal 
enrichment of the ICM , while total gas mass increases 
(|Zhang et alj|201ll ). This implies cluster star formation 
efficiency that is dependent on the halo mass, such that 
high mass halos turn a smaller fraction of their baryonic 
mass into stars. This could be the result of a constant de- 
pendence of the quenching efficiency of clusters on mass 
over all time, so that while accreting more mass from the 
field, high mass clusters preferentially return the gas to 
the IGM rather than process it into stars. Additionally, 
this difference could reflect a change in environmental 
quenching mechanisms with time. At a given epoch clus- 
ters with different halo masses have different formation 
histories, assembly times and galaxy population ages. N- 



body simulations, local stellar population studies, and 
high-redshift observations of clusters indicate that high 
mass clusters formed a larg er fraction of their stars at 
earlier times (Thomas et al. 2005) than lower mass sys- 
tems, thereby imprinting temporal variations on global 
star formation efficiencies on the cluster galaxy popula- 
tions. 

Over the mass range probed here (~10 14_15 Mq) local 
clusters show a M*/Mh a io dispersion of approximately an 
order of magnitude. The simplest interpretation of this 
is that their star formation efficiency follows the same 
ratio (with the higher mass clusters an order of magni- 
tude less efficient at forming stars than the lower mass 
systems) , at least averaged over all time. Our mea- 
sured relation of the ESFR/M c i ust0 r is consistent with 
this at ~ltr, and therefore does not require a change in 
the quenching efficiency of star formation with redshift; 
that is, the difference between high and low mass clus- 
ters can be accounted for by the relation we see below 
Z ~ 1, and does not require differences to be set in place 
at early times. According to these measurements, high 
and low mass clusters have assembled roughly the same 
amount of stellar mass since z ~ 1, though this represents 
a smaller fraction of the total stellar mass of high mass 
clusters. Thus, we confirm the fraction of a cluster's stel- 
lar mass that was formed at high redshifts increases with 
cluster mass. This is also in line with the richness depen- 
dent downsizing e ffect observe d by Gilba nTTet al.l (2008) 
(but (c.f. IPe Lucia et al.ll2008 ; Bildfc ll et alJ2012fi : lower 
mass clusters show a larger deficit of faint red sequence 
galaxies than high mass clusters to z~ 1. 

The above interpretation is in line with the proper- 
ties of low redshift halos, but it is not the only possi- 
ble explanation for Figure [7] Disentangling the effects 
of environment and galaxy mass is notoriously difficult 
and could be important here. The specific star forma- 
tion rate of individual galaxies is known to decrease with 
increasing mass even to z ~ 1 (El baz et all 120071 ) and 
thus if the mass function of in-falling galaxies varies with 
parent halo mass, such that high mass clusters preferen- 
tially accrete high mass galaxies, this would be reflected 
in the overall specific star formation rate of the clus- 
ters. This scenario is not supported by the recent mea- 
surements of the mass function of gala xies surrounding 
clusters to z ~ 1 of lVulcani et all (|2012l ) , and would nev- 
ertheless be too weak an effect to produce what is seen 
here. Still, N-body simulations do highlight a significant 
difference in the ac cretion history of halos with mass. 
iMcGee et all ([2009D show that high mass clusters form 
preferentially from group accretion, or the accretion of 
halos with M>10 13 M Q , whereas lower mass systems pri- 
marily accrete isolated field galaxies. In this scenario star 
formation is first s uppressed by the local environment 
(Balogh et al.H20iH) before accretion onto a massive halo 
and the apparent dependence on the global environment 
is simply a reflection of the mass function of accreted 
halos and the larger number of groups falling into the 
highest mass clusters. 

Finally, we recall that the RCS clusters are optically 
selected with the mass estimated through optical richness 
measurements. This means that our sample is fundamen- 
tally a stellar-mass limited sample, rather than a halo- 
mass limited sample. This could introduce a systematic 
bias in stellar-to-total mass ratio with total mass. At 
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the high mass end we will be complete for all possible ra- 
tios, but at low mass we could be biased to clusters with 
high stellar mass content, and therefore those systems 
which have experienced more efficient star formation for 
the same dark matter halo mass. This would take a flat 
SSFR/M c i us tcr relation and tilt it negatively by biasing 
the low mass clusters upwards. Such an effect should be 
recoverable provided one has an independent estimate of 
the total mass, which we currently do not. 

5.4. The radial dependence of the star formation beyond 

z> 0.5 

Figure [8] confirms the decrease in star-formation rate 
per unity galaxy toward high density regions - the 
so-called star-formation-rate density relation which has 
been well-established by many authors. While field 
studies show that the relation begins to reverse be- 
yond z ~ 1 (jElbaz et al.ll2007t iCooper et al.ll2008D . ev- 
iden ce of it reaching cluster core densities below z ~ 
1.6 (IKodama et al.|[20&l iTran et aTll2010t IMvIzzin et al.l 
I2012T) is not vet conclusive. As with the evolutionary ef- 
fect discussed above, it is exceedingly difficult to directly 
compare the results of different studies. The absolute 
value of the fraction of star forming galaxies depends 
on the depth and method of the star formation rate esti- 
mates and the method of assessing the underlying cluster 
population. Nevertheless, most studies see a smooth de- 
cline in the fraction of star forming galaxies of about an 
order of magnitude fr om the outer regions of the cluster 
(~ 9-2QO ) to the core (IKodama et ail 2004; M uzzin et al.l 
2012: Patcl et al. 2011). It is now clear that some, but 
not all, of this trend is due to the underlying mass-bias 
(|Muzzin et al.1 120121 ) ; the average mass of galaxies in- 
creases towards the cluster core and the lower sSFR of 
higher mass galaxies drives down the average star forma- 
tion rate. 

The statistics of our study are already too poor to do 
much beyond describe the overall decline (see caption of 
Figure [8]) in star formation in the cluster centers. How- 
ever, it would be interesting in future work to investigate 
the rate of decline for galaxy populations of different star 
formation rates or specific star formation rates as this 
would provide further information on the galaxy depen- 
dent properties of the quenching mechanisms. 

6. CONCLUSIONS AND FINAL REMARKS 

The simple statistical exercise presented here leads to 
a number of important conclusions: 

1. We see an steep increase in the number of IR- 
luminous galaxies (Lir > 2xlO n L0) per cluster 
mass (and by inference £SFR/M c i ustC r) of (1 + 
z )5.4±i.9 over the range 0.3 < z < 1.0. This 
evolution is in agreement with that estimated by 
other IR studies of gala xy clusters ()Haines et al.l 
120091 : IPopesso et aLll2012[ ) and shows the same level 
of increase as seen in extinction c orrected opti- 
cal s tudies to similar SFR depths (Muzzin et al. 
I2012T ). Assuming the optical population is con- 
tained within the IR population (this has not been 
demonstrated), this indicates that the optical stud- 
ies are not missing a significant fraction of dust en- 
shrouded activity. 



2. We show that the above rapid evolution can be 
accounted for entirely by the evolution in the in 
falling field population. Moreover, we show that 
the amount of IR-traced star formation per unit 
halo mass in clusters and the field are consistent; 
once normalized for halo mass clusters do not show 
an excess or in star formation deficit relative to the 
field. 

3. The IR luminous in-falling population seen here 
to z ^1 can account for 5-10% of the total stel- 
lar mass in massive clusters today. Until we have 
an estimate of the stellar mass function of these 
galaxies, it is not certain how much of this popu- 
lation contributes to the build-up of the faint-end 
of the red-sequence observed in clusters over this 
timeframe. 

4. Averaged over all redshifts, we see a decrease in the 
SSFR/M c i us tor with increasing galaxy richness of 



ESFR/M, 



lustc 



;M; 



-1.5±0.4 

lustc 



This means that the 



SFR of an individual halo decreases more sharply 
with time than the simple number count analysis 
implies. The relation is consistent with a constant 
dependence of quenching efficiency (but docs not 
require it) on halo mass over all time and can re- 
produce the dependence of the stcllar-to-total mass 
ratio on mass seen in local clusters. 

5. The radial distribution of IR galaxies in clusters 
is flat, with perhaps a slight increase towards the 
centers. This is seemingly in agreement with other 
studies which claim enhanced star formation in the 
centers of high redshift clusters. This, however, is 
not a correct interpretation; once the underlying 
galaxy density is taken into account (a step often 
neglected by other studies) we see a decrease in the 
number of IR-luminous galaxies toward the cluster 
core for z > 0.5 clusters. Thus, the SFR-density 
relation persists to cluster core densities to z~ 1 . 

As noted in the introduction the broad statistical tech- 
nique limits the conclusions to general statements con- 
cerning the average cluster and galaxy populations, and 
cannot provide information about the properties of the 
individual galaxies. Moreover, we cannot control for the 
important interdependencies of stellar mass and environ- 
ment. Such work will be the focus of later papers by 
our team and will better elucidate the history, charac- 
teristics, and fate of galaxies recently accreted into the 
cluster environment. 
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APPENDIX 

POSSIBLE SOURCES OF SYSTEMATIC ERROR IN THE NUMBER COUNTS 

In interpreting the apparent increase in IR galaxies in clusters with redshift, one must be wary of other systematic 
trends with redshift (inherent or introduced) which might mimic such evolution. There are three primary sources 
of possible systematic error in this analysis: (i) the background subtraction (ii) the use of a single template SED in 
inferring the observed 24/im flux for a given IR luminosity limit and (iii) the calculation of the r2oo radius. We discuss 
each in turn, but conclude that none are sufficient to account for the trend we see. 

Background Subtraction: Our analysis method is identical in the cluster and COSMOS fields and thus a redshift- 
dependent offset in the counts would require a systematic difference between the cluster and SCOSMOS fields that 
varies with source brightness. A difference in completeness depths could mimic source density evolution; however, a 
number density excess at high redshift would require the SCOSMOS imaging to be shallower than the cluster fields 
and we ha ve shown in Figure [3] that this is not the case; the adequate depth of the COSMOS data is further confirmed 
elsewhere (jSanders et alj|2007n . 

If, however, the SCOSMOS G03 field is located in a region of the sky of true galaxy under density at all relevant 
luminosities, this could introduce a trend similar to that observed. Again, our analysis of the number counts using our 
own photometry indicates that the number density of 24/im sources in SCOSMOS is in good agreem ent with other 
fields, and this is confirmed by the very careful analysis of the SCOSMOS team (|Le Floc'h et a l. 2009). 

A final possibility is an increased level of asteroid contamination in the RCS fields. We have not removed 
asteroids in either the RCS or SCOSMOS images, thus if the RCS-MIPS fields contain a higher number of asteroids 
than SCOSMOS, with the number density of asteroids increasing with decreasing flux, this might mimic redshift 
evolution. We have checked this by plotting the excess IR counts against the ecliptic latitude of each cluster 
and see no correlation. Moreover, the SCOSMOS field lies at lower ecliptic latitude than any of our clusters and 
therefore should have more asteroids. Therefore, it does not appear that our analysis method or inherent differ- 
ences between the cluster and field catalogs are likely to create a false increase in IR galaxies in the higher redshift bins. 

Luminosity Calculations: We adopt the Chary & Elbaz (2001) and Dale & Hclou (2002) libraries of IR luminous 
galaxy templates to predict an observed 24/im flux for Lir > 2x 10 11 L q for each cluster and are therefore applying 
the same SED to all redshifts. The Chary & Elbaz SED library has been demonstrated to accurately reproduce the 
total IR luminosity of field galaxies below z < 1.4, particularly for galaxies with L < 10 12 L Q which dominate our 
number counts, even when extrapolated from a single 24/im measurement (Magnelli et al. 2009; Murphy et al. 2009; 
Chary 2010). Thus, field studies show no evidence of a strong evolution in the IR SED over this redshift range and we 
therefore do not expect this to be a substantial effect. Nevertheless, if the average SED of IR galaxies within clusters 
changes with redshift we could be measuring the number counts to different intrinsic luminosities at different cosmic 
epochs. To reproduce the observed trend of increasing numbers of IR galaxies at high-redshift requires an enhancement 
of 24/im flux for a given Lir in the higher redshift bins. An obvious means of accomplishing this is through increased 
AGN contamination as the presence of hot dust surrounding an AGN would produce excess MIR emission thereby 
boosting galaxies above the luminosity limit. 

To investigate this possibility we have computed the IRAC MIR colours of the 24/xm galaxies in each cluster field 
to the correspon ding luminosity limit); the MIR has been used to roughly classify galaxies as possible dusty AGN 
Lacv et all 120041 ) through their location on an IRAC color-color plot. We measured the fraction of 24/zm galaxies in 
the cluster fields which have colors consistent with AGN according to the Lacy et al. classification in three separate 
redshift bins above z = 0.5 and find no change in this fraction with redshift. However, given the small number of 
excess IR galaxies in each cluster field (2-20 galaxies per cluster), relative to the high source density of the background, 
color differences will be difficult to detect using this statistical method and the results are therefore not strongly 
conclusive. This will be investigated in later work using more extensive imaging and spectroscopy. We further note 
that a change in the IR SED with redshift, due to increasing AGN or other physics, simply means the evolution in 
number counts cannot be interpreted as an increase in the total SFR; it is still a real and interesting evolutionary trend. 

Adopted Radii: Since we are counting galaxies within a given radius for each cluster and the density of IR galaxies 
themselves may not be constant with cluster-centric distance, the accuracy of the analysis also relies on the adoption 
of a consistent radial cut, taken here to be r2oo, or the approximate virial radius. We calculat e r^vin from the observed 
galaxy richness N ro d, which is an indicator of mass reliable to ~30-50% (|Rozo et al. 2009; W hite et al.l 2010, Gilbank 
et al. in preparation). We choose not to attempt to correct for incompleteness in the highest redshift bins by assuming 
a constant L F, since we know tha t the faint end of the red-sequence becomes increasingly depopulated towards higher 
redshift (e.g. iGilbank et al.l 120081) in a way which approximately mimics incompleteness. 

Systematic offsets in the measured N ga i with redshift could introduce a corresponding trend in the number counts. 
For example, if the galaxy population in the cores of clusters evolves substantially between 0.2 < z < 1.0 an optically 



14 



Webb et al. 



estimated N ro( j is unlikely to have a uniform conversion to total cluster mass. In this case we would expect the same 
N rc d to correspond to a more massive system at high redshift (because the galaxies are bluer), compared to low redshift. 
This would lead to an underestimate of r2oo and therefore a decrease in the corresponding number counts, but also an 
underestimate of the mass normalization, thereby boosting the measurement. However, initial spectroscopic results 
(Gilbank et al. in preparation) do not show a significant effect. 

The counts and mass normalization depend on different powers of the richness measure, N rc( j. As outlined in §3, 
we determine T2oo by first transforming N rc d to a velocity dispersion (a), following N Ic d ex a ; the radius is then 
linearly proportional to a. Assuming a flat radial distribution of IR galaxies (§4.5) the absolute number counts will 
increase with the area, or as r^o- The mass is a slightly stronger function of radius, rising as rfoo- These competing 
dependencies result in a bias upwards in the number counts per unit mass, but the effect cannot reproduce the evolution 
for any reasonable error in richness. For example, a systematic underestimate of the optical richness in the z — 0.9 bin 
of ~30% would require a correction of the counts downward by ^10%, preserving a strong evolutionary trend. Thus 
a large systematic in richness measurements with redshifts can steepen the evolution, but cannot completely account 
for it. 

We have further checked this effect using the IRAC imaging for the z > 0.5 sample. We carry out a similar number 
count analysis as in the rest of the paper where we measure the number of galaxies within a given radius and above 
a uniform luminosity limit, and perform background (line-of-sight) subtraction. In this case we have chosen 0.5 x 
r2oo, to maximize S/N, and Mr > M^- + 1, to remain well above the IRAC limiting depths at all redshifts. This 
is not a perfect reproduction of the MIPS analysis method since the optical/NIR galaxies are unlikely to have the 
same radial distribution as the IR galaxies (§4.4) at 0.5 r2oo, and the luminosity conversion method is different, but 
it should provide a handle on the importance of large systematic errors in the counts due to systematics in the radii. 
Within small radii, IRAC traces the stellar mass of the early-type galaxies dominating the core and should show very 
little evolution with redshift (in particular since these systems have been selected by the existence of an old optical 
population). Thus, any large evolutionary trends in the IRAC counts may be attributed to systematic effects. In fact, 
figure |9] shows no trend of increasing optical galaxy counts with redshift, in stark contrast to the MIPS analysis. 

RAW COUNTS 

Here we show the background-subtracted 24/^m counts in constant radii of IMpc and without normalizing for the 
inferred cluster mass. This plot makes no assumptions or corrections for the differences in counts which are due 
to differences in the cluster mass and does not infer an r2oo radius. It should therefore be free from uncertainties 
introduced through errors in the assumed optical-richness-mass relation, or its intrinsic scatter. On the other hand, 
working at a fixed radius and not accounting for mass should result in a larger scatter - assuming the corrections 
are properly applied. This is because IMpc corresponds to a larger fraction of r2oo for poor clusters, compared to 
richer clusters, and because higher richness clusters may have larger numbers of all galaxy populations - including the 
IR-luminous systems. The increased number of IR-luminous galaxies in higher mass clusters appears robust in that it 
is present in both the uncorrected raw data, and (at a higher level) for the corrected data. For the mass-uncorrected 
data shown here the evolution goes as (1 + z) n where n =4.4±1.9, compared to n =5.1±1.9 for the corrected version. 
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Fig. 9. — The number of IRAC detected galaxies (3.6/jm or 4.5/i4m) within 0.5 T2oo an d brighter than M*+l, after background subtraction. 
The black points correspond to individual clusters and the yellow points show the weighted average within three redshift bins. 
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Fig. 10. — The number of 24/xm galaxies above an inferred luminosity of 2X20 11 Lq within the same radius of IMpc for all clusters. This 
plot is the same as the upper panel in Figure 4, however it does not include the mass normalization and does not alter the counting radius 
according to the cluster richness. The solid line is the best-fit power law function of N ~ (l+z) 44±1,9 . 



